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Abstract 



The talk briefly reviews the current state of the art in doing multiloop QCD 
calculations in a completely analytic way. In particular, we discuss recent 
analytical calculations of 4-loop 0{o? s ) gluino contribution to R(s) and its 
implications to T tot (Z — > hadrons) and r(r~ — > v T + hadrons); 
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1 Introduction 



The firm knowledge of radiative corrections (r.c.) is a must in confronting the Standard 
Model with experimental data. While the electroweak r.c. are usually required to be 
known in the leading and next-to-leading approximations, this is certainly not always the 
case for the strong interaction ones. The current (relatively large) value of the strong cou- 
pling constant a s [Mz) = 0.1202 ± 0.0033 fl| and the remarkably high precision achieved 
in high energy e + e~ experiments make the latter sensitive to really higher order QCD r.c. 
such as 0(a^) contributions to R(s) and T tot (Z — > hadrons). 

The talk briefly reviews the current state of the art in doing multiloop QCD calcu- 
lations in a completely analytic way. It is not intended as a sort of exhaustive review of 
the vast subject of the QCD r.c. We deal with analytical calculations of 3- and 4-loop 
r.c. to 2-point correlators of bilinear quark currents. Even more, only massless correlators 
are considered (high energy limit). The choice is dictated in part by my personal taste 
and in part by the the fact that most advanced results (at least as for the number of 
loops involved) have been obtained in this particular field. The reader who is looking for 
a broader exposition including leading and next-to- leading r.c. should consult numerous 
reviews 

2 Feynman Integrals up to 3 loops 

Here we will discuss briefly the tools now available to analytically compute massless prop- 
agators in higher orders. We limit ourselves to these rather restricted class of Feynman 
integrals due to the following reasons: 

1. Practice shows that in many cases the methods of asymptotic expansions of Feyn- 
man integrals (for a recent review see ||) do produce results numerically very well 
approximating the exact results when the latter are available. These methods reduce 
initial multi-scale Feynman amplitudes to combinations of massless propagators and 
massive tadpoles. 

2. A number of problems can be eventually reduced to evaluation of massless prop- 
agators. A important example is the evaluation of so-called RG-functions (that is 
beta-functions and anomalous dimensions). 

3. Thanks to the intrinsic simplicity of the integrals under discussion, — they depend 
on only one nontrivial scale: an external momentum — their analytical evalua- 
tion proves to be feasible in quite high orders in the coupling constant. The same 
simplicity provides the possibility of constructing regular algorithms for evaluating 
these integrals as well as dedicated computer programs allowing to perform the 
calculations in a convenient and automatic way. 

For brevity massless Feynman integrals depending on exactly one external momentum 
will be denoted by p-integrals. At the moment there are tools to analytically compute 
arbitrary one- two- and three-loop p-integrals (see below). Fortunately, in many impor- 
tant cases one is interested only in the absorptive part of massless two-point correlators. 
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Figure 1: Some p-integrals: (a) the generic one- loop p-integral and (b) an example of 
primitive five-loop p-integral; (c) the master two-loop p-integral. 



In this case available theoretical tools are enough to guarantee at least in principle the 
analytical calculability of absorptive part of an arbitrary 4-loop p-integral (see below). 

One-loop p-integrals 

We start from a well-known elementary formula for a generic 1-loop p-integral (see 
Fig. [l]a; we shall consider Feynman integrals in the Euclidean momentum space throughout 
this section) 
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It is of importance to note that any p-integral depends homogeneously on its external 
momentum. This facts allows the immediate analytic evaluation of the whole class of 
primitive p-integrals which, by definition, may be performed by repeated application of 
the one-loop integration formula. For example, the five-loop scalar integral of Fig. |l]b is 
performed by (P with the result 



(Vr(47r 2 ) 2 -)) 5 (g 2 r(G(l, 1)G(1, e)) 2 G(l + 2e, 1 + 2e) . 



(2) 



Two-loop p-integrals 

Not all p-integrals are primitive ones. One first encounters nontrivial p-integrals al- 
ready at the two loop level. While one-loop integrals are performed with ease the eval- 
uation of the master two-loop diagram (see Fig. [Tjc) is not trivial. The corresponding 
Feynman integral reads 



(47r) 4 - 2t F( ai 
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with the loop momenta 
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A closed expression for the function F(a\, . . . a§) for generic values of the arguments is 
not known. However, results do exist for particular cases. The first one, valid for a 
generic space-time dimension D, was obtained with the help of the so-called Gegenbauer 
polynomial technique in x-space (GPTX) ||. It reads 

x {a[G(a + 1,(3)- G(a + 1, (3 + e)] + (a «-> (3)} 

It has been also shown in Ref. || that similar results may be obtained for the case when 
the indices a 2 , a 3 and a 5 are integers while cti and a 4 are arbitrary. 

In practice one often needs only a few first terms of the expansion of F(ai . . . a^) in 
the Laurent series in e. This expansion is known for generic values of the ot\ . . . ot§ up to 
a fixed (quite high) order (see Refs. || |10[ and references therein) . 



Three-loop p-integrals 

In principle GTPX is also applicable to compute some non-trivial three-loop p-integrals. 
For example, the basic scalar non-planar three-loop diagram of Fig. |3|a was first calculated 
via GPTX in Ref. ||. However, calculations quickly get clumsy, especially for diagrams 
with numerators. 

The main breakthrough at the three-loop level happened with elaborating the method 
of integration by parts of dimensionally regularized integrals Refs. |7], |J. The key identity 
for the method is£] 

/d°<A I( *,... )30 , (5) 

where I(£, . . .) is a Feynman integrand and £ is one of its loop momenta. The identity 
reflects the possibility of neglecting the surface terms, which holds true in dimensional 
regularization The use of ([5]) along with tricks like completing momentum squares 
and cancelling similar factors in the nominator against those in the denominator consti- 
tutes the essence of the approach. The identity depicted in Fig. ^| is a typical example of 
relations obtainable with the help of the integration by parts method. 

It should be, however, heavily stressed that, the validity of such operations for divergent 
dimensionally regulated integrals with deeply intermixed UV and IR (sub) divergences is 
not obvious any more. It had to be rigorously justified within a proper generalization of 
the dimensional regularization itself. It has been done in Ref. |12[ (see also a monograph 
II)- 

The general scheme of the use of the integration by parts method is based on the 
exploitation the identities of type (fD in the form of recurrence relations with the aim to 
express a complicated diagram through the simpler ones. All (about a dozen) topologically 
different three-loop p-integrals were neatly analyzed in Ref. || and a concrete calculational 
algorithm was suggested for every topology. As a result the algorithm of integration by 
parts for three-loop p-integrals was developed. The algorithm constitutes a series of 
involved identities which are used to identically transform any three-loop p-integral into 



x For two-loop massive integrals a similar identity was used in the classical work by 't Hooft and 
Veltman Ref. El. 
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«<D = 0-oc> 



Figure 2: The exact relation expressing a nonprimitive two-loop scalar p-integral through 
primitive integrals; a dot on a line means a squared scalar propagator. 




a sum of primitive one-loop p-integrals and two basic three-loop p-integrals pictured in 
Fig- § 

It goes without saying that the calculation of higher order corrections in gauge theories 
is almost impossible without intensive use of computer algebra methods. In addition to 
the old problem of taking long traces of Dirac 7 matrixes, the algorithm of integration 
by parts, when applied even to a single 3-loop p-integral, generically produces dozens or 
even hundreds of terms. At the moment there exist essentially three different packages 
which implement the algorithm. For p-integrals they are written in SCHOONSCHIP 
(see Refs. [[15], @ and in FORM g7| (see Refs. @, 0). 



3 Infrared Rearrangement and i?*-operation 

At the moment there is no any general algorithm allowing to analytically compute arbi- 
trary 4-loop p-integrals. The problem of creating such an algorithm seems to be hopelessly 



difficult. See Ref. [2(J where the point is dealt with in some detail and also a discussion 
in the conclusion of the present talk. 

Nevertheless, a large group of important 4-loop problems — calculation of RG-functions 
or, equivalently, UV renormalization constants — proves to be reducible to the 3-loop case 
and, thus, be doable in a completely analytic way. This is because in MS scheme any UV 



counterterm is polynomial in momenta and masses [[H]]. This observation was effectively 



employed in Ref. p3] to simplify considerably the calculation of UV counterterms. The 



method was further developed and named Infrared Rearrangement (IRR) in Ref. Q. It 
essentially amounts to an appropriate transformation of the IR structure of FI's by set- 
ting zero some external momenta and masses (in some cases after some differentiation is 
performed with respect to the latter). As a result the calculation of UV counterterms is 
much simplified by reducing the problem to evaluating massless p-integrals. The method 
of IRR was ultimately refined and freed from unessential qualifications in Ref. |23| . The 
following Theorem has been proven there by the explicit construction of the corresponding 
algorithm: 

Any UV counterterm for any (h+l)-loop Feynman integral can be expressed 
in terms of pole and finite parts of some appropriately constructed (h)-loop 
p-integrals. 

In many important cases one is interested only in the absorptive part of massless 
two-point correlators. In this case available theoretical tools are enough to guarantee 
at least in principle the analytical calculability of absorptive part of an arbitrary 4-loop 
p-integral. This is eventually due to the fact that the absorptive part of a four-loop p- 
integral is fully expressible in terms the corresponding four-loop UV counterterm along 
with some three-loop p-integrals. As a typical example we will consider in next section 
the ratio R(s) = o"tot( e+e ~ — > hadrons)/a(e + e _ — > which is essentially given by 

the absorptive part of the vacuum polarization. 

The Theorem coupled with the the integration by parts method solves at least in 
principle the task of analytical evaluation of RG functions and absorptive parts of 4-loop 
p-integrals. It should be noted that the -Reoperation is essential to prove the Theorem, 
though in most (but not in all) practical cases one could proceed without it. However, 
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such a practice, in fact, forces diagram-wise renormalization mode, what, in turns, brings 
down a heavy penalty of manual treatment of hundreds of diagrams. 

Indeed, in genuine four-loop calculations the reduction to three-loop p-integrals is far 
from being trivial and includes a lot of manipulations. Typical steps here are 

a to reduce the initial Feynman integral to logarithmically divergent ones via a proper 
differentiation with respect masses and external momenta; 

b to identify UV and IR divergent subgraphs of the resulting integral; 

c to remove in a recursive way the corresponding UV and IR divergences; 

d to compute resulting p-integrals. 

Among these steps only the calculation of p-integrals can be at present completely per- 
formed by a computer. All others, especially b and c are difficult to computerize. As 
a result, in spite of the fact that the necessary theoretical tools have been around for 
more than a decade, yet until very recently there existed just one QCD four-loop result: 
the 0(a? s ) contribution to R(s). It was done in Refs. [24, 25] within the diagram-wise 
renormalization approach. The core of the problem is the calculation of the four-loop 
contribution to the photon anomalous dimension entering into the RG equation for the 
photon polarization operator. The initial 98 four-loop diagrams contributing to the photon 
polarization operator proliferate to about 250 after the IRR procedure is applied. In ad- 
dition these diagrams contain about 600 various subdiagrams which should be computed 
separately in order to subtract UV sub divergences. No wonder that both calculations 
were done in a particular gauge — the Feynman one — and, thus, a valuable possibility 
of using the gauge invariance as a strong test of the correctness of the results had not 
been used at all. 

The calculation of Ref. did not use the .Reoperation at all while that of Ref. [E3] 



employed it only for a few diagrams. We will see in the next section that a proper use of 
power of the -Reoperation allows one not only to repeat with ease these calculations in 
the general covariant gauge but also to extend them considerably. 



4 4-loop case: photon vacuum polarization 

The photon vacuum polarization function IT(— g 2 )is defined through the correlator of the 
hadronic EM current = J2f Q pl> fl^f as follows 

n w (?) = (4vr) 2 * J dxe^(0\T[ j^(x) 3 T(0) P = -3g 2 n(-g 2 ) . (6) 

The optical theorem relates the inclusive cross-section and thus the function R(s) to the 
discontinuity of II in the complex plane 

R(s) = — hnU(-s-iS). (7) 

4:71 

The renormalization mode of the polarization operator II (Q 2 ) reads 

n(Q7/Aa s ) = Z em + n (Q 2 ,a°),, (8) 



6 



where a s = g 2 /(47r) is the strong coupling constant. Within the MS scheme (here and 
below we are using a convenient combination a s = a s /7i) 

z cm = E (z™)^, (9) 

with the coefficients (Z em )^. being pure numbers and D = 4 — 2e standing for the space- 
time dimension. 

As a result we arrive at the following renormalization group (RG) equation for the 

2 

polarization operator (Lq = In ^ ) 

^-U = Q 2 7 em(a s ) - (^{a s )a s ^j II. (10) 

Here the photon anomalous dimension and the /3(a s )-function are defined in the usual 
way 

7em = [i 2 -f^(Z cm ) - eZ cm = ~53(z + l)(Z cm )na l s , (11) 
^A_ as = as p( as ) = _ £ ^\ (i 2 ) 



The relation ( ]T0| ) explicitly demonstrates the main computational advantage of finding 
first the polarization function H(Q 2 ) against a direct calculation of R(s) in the case of 
massless QCD Indeed, in order a™ the derivative gf^n and, consequently, R(s) depends 
on the very function II which is multiplied by at least one factor of a s . This means that 
one needs to know II up to order a™ -1 only to unambiguously reconstruct all Q-dependent 
terms in II to order a", provided, of course, the beta function and anomalous dimension 
7e m is known to order a". On the other hand, as we have discussed before the calculation 
of an anomalous dimension or a beta-function is known to be much easier than computing 
a correlator of the same order in the coupling constant. 

Thus, in order to check the results of Refs. [24], |25| by an independent calculation one 
should compute II in order a 2 s (a dozen three-loop diagrams) and 7 em in order (more 
than a hundred of four-loop diagrams). The first task is almost trivial nowadays (see 
the previous two sections) so let us concentrate on the second one. As we have already 
discussed, the problem is doable via IRR and the integration by parts method but the 
amount of various non-trivial manipulations with separate diagrams to be performed by 
hand is really terrific. 

In a recent work |2(J it has been demonstrated how the formalism of the .Reoperation 
can be applied to derive an explicit formula for the renormalization constant Z em and, 
consequently for 7 em . The formula reads: 

Z em =-K e ^^Tr[6Tl(0,0,a° s ) la }Z^- 5 -^Z^ (13) 
- 2D{ 1 _ 1} □ g Tr[ 7a G°(p + q, ag)r° (p, q, a° s )G°(p, a° s )]\ q = Q } ■ 
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Here G° and T° are the full quark propagator and the EM current vertex function re- 
spectively; the integration with respect to the loop momentum p with the weight function 
iw is not explicitly displayed. Z 2 is the quark wave function renormalization constant. 
Several extra comments are in order regarding this formula. First, Eq. ( |13"1) is, rigor- 



ously speaking, applicable as it stands only to the so-called non-singlet diagrams, that 
is to those where both EM currents belong to one and the same quark loop. The four 
singlet diagrams, violating this requirement, appear first in order and should be treated 
separately. 

Second, by ST^(p, q, a°) we denote the vertex function of the electromagnetic current 
with the tree contribution removed. The "tilde" atop the index a again that in every 
diagram the quark propagator entering to the vertex j a is softened at small momenta by 
means of the auxiliary mass mo according to the rule 

la -> la = la V % I (p 2 ~ TT%). (14) 

The bare coupling constant a° is to be understood as a s = Z a a s , with Z a being the 
coupling constant renormalization constant. 

An inspection of (p~3|) immediately shows that, in order to find the (n + l)-loop correc- 
tion to Z em , one needs only to know the renormalization constants Z 2 and Z em to order 
a™ as well as the bare Green functions 

G°( P ,a% ^[r£(p,ff,a?)]| g = , □ 9 [I>,g,a°)]| (? = , 5r°(0,0,a°) (15) 

up to (and including) rz-loops, that is to order (a°) n . 

We have computed with the program MINCER [18[] the unrenormalized three-loop 
Green functions (|T5|) as well as the quark wave function renormalization constant Z 2 to 
order a 3 . The calculations have been performed in the general covariant gauge with the 
gluon propagator {g^ u — (,^r L )/q 2 . We have also taken into account the singlet diagrams 
as well as extra diagrams with some of virtual quarks replaced by colour octet neutral 
fermions. In the minimal supersymmetric standard model such a fermion known as gluino 



appears as the superpartner of the gluon [EH 



Our result for R(s) with /i 2 = s reads [26 



R(s) = 3^Q 2 f {l + a s + a 2 s r 2 + a 3 s r 3 } + al[J2Qf) rf, (16) 



/ 




where rf = (|§ - §C(3)) and 



365 . . 11 2 , . 11 

T2 = — -HC(3)--n / + -C(3)n / - T% + 2C(3) % , 

87029 121 2 1103 A/ , 275 A/ x 7847 11 2 

r 3 = 7r 2 C(3)H C(5) nf H n 2 n f 

3 288 48 4 sv ; 6 SV ; 216 f 36 ' 

262 25 151 a 1 2 2 19 2 

+ C 3 )n f C{h)nt-\ n f ti n f ( 3 n f 

9 SV ; t 9 SV ; ' 162 f 108 f 27 SV ; f 

32903 11 „ 277 u/ , 25 151 



— n-rig + — C(3)%- -^-((5)% + — n f n 



288 9 12 9 3 " v ' 9 3 " v ' 9 27 

1 2 38 . , 151 2 1 2 2 19 //o\ 2 

— 7T n f n 5 r(3)n/-n 5 H n- tt n- C(6)n~. 

18 ; 9 9 18 9 12 9 3 V ; 9 
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Here rig is the number of neutral gluino multiplets which we take either zero or one. We 
observe that neither 7 em no R(s) depend on the gauge fixing parameter £ as it must be. 
If rig is set to zero then R(s) is in complete agreement with the results of Refs. [24, |25[ . 
In the numerical form, 

R( s ) = 3 Q) {l + a-s + al (1-98571 - 0.115295™/ - 0.345886n 5 ) 
/ 

+a 3 s (-6.63694 - 1.20013n/ - 0.00518nJ - 2.85053n 5 
-0.03107n f n g - 0.04661n?) } - a 3 s (j^ Qfj 1-2395. 

5 0(a>l) gluino contribution to 

T tot {Z -> h's) and r(r" i/ T + h's) 

The result (|16"D can be straightforwardly applied to find the 0(a;f) gluino contribution 
to the Z-boson decay rate into hadrons (I^)- As is well-known this decay rate may be 
viewed as an incoherent sum of vector (T^) and axial (T^) contributions. For massless 
u, d, s, c and b quarks the QCD corrections to Y v z are in one-to-one correspondence to 
those for R(s). If one neglects as we do 2 the power suppressed terms of order Mf /(4M 2 ) 
and higher the same statement is valid also for T z except for a specific subset of so-called 
singlet diagrams. As a result the total decay rate is naturally presented as a sum of three 
terms 

r A,S _ riV , r A,NS | r A,S /-. 7 \ 

L z — L z L z ' 1 z ■ \ 11 ) 

Here T^' s and T^' NS stand for the contributions to due to singlet and the rest (non- 
singlet) diagrams respectively. The first two terms in (|17D are directly expressed through 
the coefficients r\ — r% and rf as follows: 

z+ z J '^" 2 ^ 1 + a s + a 2 s r 2 + a 3 s r 3 





r 

where r = G f MI/2AV2tt = 82.94 MeV, v f = 2l( - 4Q f sin 2 6 W , a f = 2l(, a s = a s (M z ) 
and we have set \x = Mz- 

The singlet diagrams are sensitive to the huge mass splitting in the top-bottom dou- 
blet and should be computed afresh. To order a 2 this was done in Ref. [Q (even without 
neglecting the power suppressed terms) while the 0(a 3 ) corrections were computed in 
|29|| . In fact, no extra calculations are necessary to get the gluino contribution. A sim- 



ple inspection of relevant diagrams immediately reveals that gluinos appear exclusively 
through the one-loop fermion correction to a gluon propagator This means that the result 
in QCD with gluinos can be obtained from the one in pure QCD by the replacement 
n f ~^ n f + 3n<j. It reads 



r A z s = r U 



h 3 In— | 

4 Mf 



2 These terms have been found to be extremely small in Ref. ¥Tt 
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+ a 3 s — — + — — + 3C(3) + — In— | + — In 2 — | (18) 




19 Ml 33 2 Mf 

2- ln W + T ln W 

+ ^ + 3 ns )|^-^ + -ln s |--l„^ 

where M t is the pole mass of the top quark and rif = 5 is the number of (light) quarks 
flavours. Numerically, the gluino term is tiny: the above equation with n/ = 5, and 
Mf/Mf = (91.187/180) 2 is evaluated to 

r^' 5 = T {-13.33a 2 + a 2 [-52.514 - 2.8197^]} . (19) 

Another obvious application of Eq. ( JT5| ) is the gluino contribution to the semihadronic 
decay rate of the r-lepton. As is well-known in the massless limit the perturbative con- 

•i ! • n T(r — > u T + hadrons) . in . 

tnbutions to the ratio R T = — can be written as follows (for a review 

r(r -> v T e u e ) 

including non-perturbative and other effects see Ref. 

r M T ds 

Rr = 2j o — (1 - ,/M 2 ) 2 (l + 2s/M*)R(s), (20) 

where M 2 is the r-lepton mass. Here R(s) is given by the expression (|T6[) for -R(s) 
with 3J]Q 2 replaced by 3(|Vy 2 + |Kts| 2 ), n f = 3 and (I]Q/) 2 set to zero. After a 
straightforward integration with respect to s in ( ^0|) one arrives at (for gluino in the octet 
representation) 



Rr = 3(\V ud \ 2 + \V US \ 2 ){1 +a s + a 2 [5.2023 - 1.13755n § ] 

} , (21) 



+ al 



26.366 - 21.0358n g + 1.42119n? 



with a s = a s (M T ). 

Note that the result ( PTj ) is applicable only if the gluino mass is significantly lower than 
that of the r-lepton. Taken at its face value it could spoil the current agreement between 
the values of Aqcd extracted from the R T and Y(Z — > hadrons). Such a conclusion seems 
to be premature as for a meaningful phenomenological discussion of the 0{a z s ) gluino 
contributions one should also take into account the running of the coupling constant in 
the next-next-to-leading order. This requires the knowledge of the gluino contribution to 
the three-loop coefficient $2 which is not yet available in the literature. 



6 Conclusion 



The skillful use of the -Reoperation leads to significant simplifications of analytical 4-loop 
calculations if we are interested in the divergent part of 4-loop diagrams. Thus, 4-loop 
calculations in QCD are gradually becoming practically feasible. 

It is only natural to ask now about next level, that is about possibility to analytically 
compute finite parts of 4-loop p-integrals and consequently (according to the Theorem 
from Section 3) divergent parts of 5-loop p-integrals. Even a superficial glance at the 
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problem shows that it hardly could be done completely "analytically". Still, there is, 
to my opinion, a chance that in a sense the problem could be solved. I mean that a 
better understanding of all kinds of identical relations connecting various p-integrals| 
could eventually result to the reduction of an arbitrary 5-loop p-integral to a combination 
of some limited number (a few dozens?) of master p-integrals. Once it is done there should 
be not very difficult to evaluate the latter analytically or numerically with sufficiently 
high accuracy. A fresh example of such approach has been recently demonstrated in a 
work by Laporta and Remiddi who achieved a computer-algebraic reduction and eventual 
analytical calculation of all 3-loop diagrams contributing to the electron's anomalous 



magnetic moment | 33|| . 



The present work has been submitted to the HEP-PH e-print archive today, that is 
29.10.1996. A related article 0] was posted on the same archive yesterday. Our results 



for the 0(a z s ) gluino contributions to R(s) and to related quantities have been confirmed 
there by an independent calculation. 
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